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Abstract 

We study the spectrum in such a VT— symmetric square well (of a diameter L < oo) 
where the "strength of the non-Hermiticity" is controlled by the two parameters, 
viz., by an imaginary coupling ig and by the distance £ < L of its onset from the 
origin. We solve this problem and confirm that the spectrum is discrete and real 
in a non-empty interval of g < go(£,L). Surprisingly, a specific distinction between 
the bound states is found in their asymptotic stability/instability with respect to 
an unlimited growth of g beyond go(£,L). In our model, all of the low-lying levels 
remain asymptotically unstable at the small f < i and finite L while only the 
stable levels survive near £ « L < oo or in the purely imaginary force limit with 
< £ < L = oo. In between these two extremes, an unusual and tunable, variable 
pattern of the interspersed "robust" and "fragile" subspectra of the real levels is 
obtained. 
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1 Introduction 



Around 1992, Daniel Bessis succeeded in attracting attention of a few people to a cer- 
tain toy Hamiltonian (with some relevance in quantum field theory) which appeared 
to produce the real and discrete spectrum of energies in spite of being manifestly non- 
Hermitian [1]. A few years later, Bender and Boettcher returned to his mind-boggling 
problem and published a numerical study [2] of the whole class of the perceivably 
more general one-dimensional Schrodinger equations 

4>(x) = Eij){x) (1) 

where, in our present perspective, the real component of the potential was assumed 
spatially symmetric while its Hermiticity-violating partner was chosen as spatially 
antisymmetric, 

VV(x)V = V(-x) = +V(x), VW(x)V = W(-x) = -W{x). 

The latter study confirmed that the similar models [exhibiting, obviously, the parity 
(V) times time-reversal (T) symmetry] may possess both the purely real and partially 
(or, perhaps, completely) complex spectra. The Bender's and Boettcher's Figure 1 
(loc. cit.) illustrated the existence of the spectrum which proved "robustly real", 
i.e., real in a wide range of parameters of their "massless" VT— symmetric model. 
In contrast, a merely slightly modified "massive" VT — symmetric model of their 
Figure 3 (loc. cit.) behaved quite differently. The values of many of its energy 
levels proved extremely sensitive to the very small variations of the parameters and, 
moreover, even the very reality of some energies proved "fragile" in the sense that 
after a very small change of a parameter of the model, certain energy pairs merged 
and disappeared forming, presumably, the complex conjugate pairs. At present, 
many more similar and more or less purely numerical examples exists (cf., e.g., the 
recent paper [3] for a sample of references). 
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+ V(x)+iW(x) 



The recent progress in our understanding of the various VT— symmetric quantum 
Hamiltonians H may be briefly summarized as an observation that their symmetry is 
important. Firstly, it was established that the time-reversal-type antilinear operator 
factor T merely mediates the Hermitian conjugation A — > [4, 5]. The role of 
parity V is more subtle and seems to offer the main mathematical key to the study 
of the VT — symmetric quantum Hamiltonians H within the so called Krein-space 
theory (cf., e.g., ref. [6] for a nice as well as concise introduction to this language). 

On the background of these mathematical observations, the formalism lost its 
originally highly enigmatic features in the context of physics. During the last two or 
three years, the use of the VT— symmetric quantum models has in fact been accepted 
as just opening new horizons within the standard Quantum Mechanics. At present, 
virtually all the people active in the field would agree that it is only necessary to 
make the resulting physical picture complete by a revitalization of its probabilistic 
contents and tractability. This is being achieved via an introduction of the "missing" 
(and, in fact, quite nontrivial) metric r\ ^ I in the Hilbert space of states [6, 7, 8]. 

The temporary doubts and puzzles related, typically, to the applicability of the 
formalism look, at least roughly, clarified. One feels urged to return to many recently 
neglected and apparently evasive and mathematically more subtle questions like the 
problems of the robustness/fragility of the individual energies or of a global typology 
of the spectra. We believe that it is time for their deeper and more technical study 
via, say, simplified and, first of all, non-numerically tractable models. A new one, 
with rather surprising properties and descriptive features of the spectrum, is to be 
proposed and analyzed in what follows. 
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1.1 Non-Hermit ian square- well- type models 

Within VT — symmetric Quantum Mechanics a one-parametric non-Hermitian square 
well (NSW) model has been described in ref. [9]. A key merit of the NSW model lies 
in a combination of its straightforward mathematical solvability with an exceptional 
transparency of its applications. In this way, the NSW model was able to offer an 
insight into the mechanism of the spontaneous VT— symmetry breaking [10]. Next, 
due to its elementary character, the NSW model has been selected by Bagchi et al 
[11] as a starting point of a systematic supersymmetric generation of solvable non- 
Hermitian Hamiltonians with VT— symmetry and real spectra. Last but not least, 
Mostafazadeh and Batal [7] choose the NSW model in their very recent illustrative 
application of the VT— symmetric Quantum Mechanics in its present, mathemati- 
cally as well as physically more or less consistent updated form (readers may consult 
some of the available reviews for more details [12]). 

In our recent paper [13] we revealed that a certain "hidden" shortcoming of the 
NSW model may be seen in its "fragility" , i.e., in an instability of all the higher energy 
levels with respect to a certain highly speculative form of a complex-coordinate 
perturbation. Although such an observation does not have any immediate impact 
on the applications of the NSW model in refs. [7, 10, 11], certain doubts survive 
concerning the possible manifestations of some more serious instabilities in some of 
the generalized, NSW-type (NSWT) models. 

For our present purposes let us vaguely characterize the latter NSWT potentials as 
piecewise constant. Then we may immediately recollect the existence of several "user- 
friendly" NSWT examples incorporating square-well models on a compact domain 
[14] or systems based on the use of point interactions [15]. Unfortunately, even within 
this class, the expectations concerning the stability of the spectrum are not always 
fulfilled. One may recollect, e.g., a spontaneous complexification of the high-lying 
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part of many NSWT spectra as detected in very early numerical studies of certain 
particular potentials in ref. [16]. The phenomenon looks puzzling and makes all the 
NSWT models worth a more detailed non-numerical study. 

1.2 The choice of a specific example 

In applied quantum mechanics the construction of the majority of phenomenological 
models relies quite heavily on the correspondence principle which tries to connect 
each quantum model with its classical predecessor. VT— symmetric Quantum Me- 
chanics offers a weakening of this connection [17]. The operator of parity V is indef- 
inite so that, as we already mentioned, the formalism requires an explicit additional 
construction of a Hamiltonian-dependent positively definite metric rj > in Hilbert 
space. Equivalently, this may be mediated by the construction of a quasi-parity Q 
[18] or charge C [19], both defined as a product vfV . In practical calculations this 
means that the metric is often being introduced in a suitably factorized form [20]. 

It is worth adding that the quasi-parity in r\ = QV is easily defined in some exactly 
solvable examples [18] while the charge in 77 = CV has immediate connotations in 
field theory [19]. In between these two extremes the authors of ref. [7] revealed that 
the application of the formalism to the particular NSW model proves facilitated by a 
perturbative connection between the NSW model and a Hermitian square well. Their 
construction of r]^ NSW ^ profited from the existence of a finite-dimensional matrix 
approximation of the non-Hermitian part of the NSW Hamiltonian. A transition 
to the extended NSWT class of models looks promising and co-motivates also our 
present project. 

Within such a framework we intend to pay attention to the family of Schrodinger 
equations (1) where the interaction is non-Hermitian but manifestly VT — symmetric. 
For the sake of definiteness we shall contemplate the less interesting real part of the 
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potential just in the most elementary infinitely deep square-well form, 



V(x) = < 



+00 


+00 



for < 



x > L 
-L < x < L 
x < -L. 



(2) 



This means that all our wave functions have to vanish at its walls, 



^(-L)=^(L) = 0. 



(3) 



By adding any imaginary interaction we break the Hermiticity of the Hamiltonian. 
By doing so in the VT— symmetric manner we preserve a chance and good hope of 
having the energies real [2]. 

For the sake of defmitness and in a way generalizing the NSW model of ref. [9] 
we shall assume that the Hermiticity-breaking term W is composed of two purely 
imaginary steps which both vanish inside a subinterval (—£, £) of the interval (— L, L), 



W(x) 





-ig 



for 



Re x > £ > , 
Re x e (-£,£), 
Re x < -£. 



(4) 



A priori, the strength of the Hermiticity-violating imaginary force may be expected 
proportional to the coupling g > and inversely proportional to £ < L. 

Our interest in the particular two-parametric model (4) results from the obvious 
need of an enhancement of flexibility of its one-parametric NSW predecessor and 
also from the lasting possibility of its rigorous mathematical description by means 
of the efficient moving-lattice method of ref. [13] (reviewed also briefly in Appendix 
A below). Among additional purposes of the study of the similar NSWT models 
one may list a search for reliable comparisons between different potentials revealing, 
hopefully, some new, unnoticed characteristic features of their spectra. One would 
like to understand, i.a., how the details of the shape of W(x) could influence the 



stability of the spectrum, or how one could control the domain of parameters where 
all the energies remain real. 

Some of the NSWT studies have been motivated by their potential capacity of 
mimicking the properties of unsolvable models and, in particular, of one of the most 
popular VT— symmetric toy interactions W(x) = ix 3 [21]. Some parallels are defi- 
nitely there since in the latter unsolvable case the spectrum was proved real, non- 
negative and discrete [22]. Of course, there are always good reasons for an intro- 
duction of more parameters in NSW. Thus, the new freedom of a weakening of the 
non-Hermiticity by the choice of I > might simulate analogies with the Bender's 
and Boettcher's generalized VT— symmetric family W(x) = — {ix) 3 ~^ characterized 
by an abrupt change of its spectral properties at /i — 1 and by the spontaneous com- 
plexification of all the sufficiently high- lying energies inside the interval // e (1,2) of 
the shape-parameter [2]. 

The possibility of the latter correspondence passes an easy test at £ = and 
L = oo when the general solutions of our Schrodinger eq. (1) are mere exponentials 
at any real g > 0. Once we demand that they vanish in infinity we have 



B + exp(—crx), a 2 = ig — E , IZea > 0, x E (0, oo), 

(5) 

B_ exp(cr'rr), a' 2 = -ig - E , Ilea' > 0, x£ (-oo,0). 



When x — > ± the coincidence of the right and left limit of ip(x) itself specifies the 
normalization, B + = while the second matching rule tp'(0 + ) = i^'(0~) implies 
that a = —cr f , i.e., equation (5) has no solutions at g > 0. It is of no avail to admit 
that Tie a — > and He a' — > and to employ the scattering boundary conditions 
since, unless g — 0, the matching-compatible states remain always incompatible with 
our differential Schrodinger equation on a half-line. 

We may conclude that both the discrete and continuous spectra are empty at 
t = for g > and L = oo. This re-confirms our above expectations since the 
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emptiness of the spectrum also characterizes the Bender's and Boettcher's toy inter- 
action W(x) = —(ix) 3 ~^ at the Herbst's extreme shape parameter /j, — 2 [23]. At 
the same time, the spectrum abruptly ceases to be empty at \i < 2 [2] as well as at 
£ > while L = oo (cf. the proof of this assertion as given in Appendix B below). 



2 The method 



2.1 Wave functions and their matching 

As long as our potential is piecewise constant atO<£<L<oowe may postulate 

4>-(x) = B_ sinh/t*(L + x), xe(-L,-£), 
i(>(x) = { xj, Q {x) = C coskx + iD sinkx, x e (-£,£), ( 6 ) 
ip+(x) = B + sinh k(L — x), x e (£, L) 

where k = s+it, E = k 2 = t 2 — s 2 , g = 2st > and where s, t and k are assumed real 
and, for the sake of defmiteness, positive. In the other words, we assume that within 
a not yet specified non-empty domain of parameters g and £ the VI symmetry of 
the wave functions remains unbroken. In the way proposed in ref. [9] we prescribe 
the phase, 

ip(x) = real symmetric + imaginary antisymmetric 

and deduce that C and D are real. Next, we differentiate 

i//_(x) = k*B_ cosh/t*(L + x), x e (-L, -£), 

ip' (x) = —k C sin k x + i k D cos kx, x e (—£, £), 
ip' + (x) = —k B + cosh k (L — x) , x G (£, L) , 



ip'(x) 



7 



and write down the following four matching conditions, 

=-0 o (-f), i.e., B_ sinh/t*(L - £) = C cosk£-iD sink£, 

ij/_(-e)=il/ (-e), i.e., k*B„ coshn*(L-£) = kC sinki + ikD cosk£, 

ip + (£) = ip (£), i.e., B + sinh k (L - £) = C cos k £ + i D sin k £, 

ip+(£) = 4>o(£), i.e., -kB + cosh k(L-£) = -kC sink £ + ikD cosk£. 

Two of them define the (complex) values of B± so that we are left with the pair of 
the matching constraints, 

(kC sink£ + ikD cos k£) sinh.K*(L-£) = (C cos k£ - i D sink£) k* cos\ik*(L-£) 

(kC sin k £ — i k D cos k £) sinh k(L — £) — (C cos k £ + i D sin k €) k cosh k(L — £) . 

These two relations are complex conjugate of each other so that we have to consider 
just one of them, say, 

(kC sink£-ikD cosk£) sinh(s + it) (L — £) — 



with k > 0. 

2.2 Matching equations in the a — r — g space 

After we abbreviate a = s (L — £), r = t (L — £) and g = k£, equation (7) reads 

g(L — £) (C sin g — i D cos g) [sinh a cos r + i cosh a sin r] = 



(C cos k £ + i D sin k£) k cosh(s + it) (L - £) . 



(7) 



£{a + ir) (C cos g + i D sing) [cosh a cos r + i sinh a sin r] . 



(8) 



We have to keep in mind that 
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while the respective real and imaginary parts of eq. (8) have to be treated as inde- 
pendent equations 

g(L — i) (C sin g sinh a cos r + D cos g cosh a sin r) = 

— £ [a (C cos g cosh a cos r — D sin g sinh a sin r) — 
— r (C cos £> sinh <x sin r + .D sin £> cosh c cos r)] (9) 

and 

g(L — £) (C sin £> cosh a sin r — D cos £> sinh a cos r) = 

= £ [a- (C cos g sinh a sin r + sin £> cosh a cos r) + 

+ t (C cos £ cosh a cos r — D sin £ sinh a sin r)] . (10) 

In the next step we notice that the latter equations form a linear algebraic homo- 
geneous set for the two coefficients C and D. They possess a nontrivial solution if 
and only if the secular determinant V vanishes. After we abbreviate Q = tan g (= 
a quickly oscillating function of g), T — tanr (= a quickly oscillating function of 
r) and £ = tanhcr (= a monotonous and bounded function of a) we can evaluate 
T>. After a lengthy calculation the secular condition T> = acquires the following 
compact form 

X(a) + Y(t) + F(R) [x(a) + y(r)] = (11) 



where 



1 + S 2 

X(a) = — - a 2 = a 2 cosh 2 a , 

1 - T 2 

y (r) = t 2 = t 2 COS 2 T , 

v ; 1 + T 2 

x ( a ) = i_^2 a= 2 fTS1Ilh2ff ' 

T 1 

v(t ) = — r = - r sin 2 r , 

y ' 1 + T 2 2 



We may re-scale our coupling g = 2 Zj (L— £) 2 and conclude that our Z— independent 
secular equation (11), 

sin 2 g(R) \a 2 cosh 2 a + r 2 cos2r] + i?cos2 g(R) [a sinh2o- + r sin2r] = (12) 

only has to be complemented by the two trivial constraints 

ar = Z, T 2 -a 2 = R 2 . (13) 

The triplets of roots R n , a n and r n of this triplet of equations with n — 0, 1, . . . define 
all the bound-state energies E n by the elementary formula 

En = (L - £) 2 Rn = (L - £) 2 _ ^"^ ' 
In an indirect check of the recipe we may recollect its £ — > (i.e., g — * 0) limit and 
conclude that our present eq. (12) degenerates smoothly and correctly back to the 
known secular £ = equation {cf. eq. Nr. (9) in ref. [9]}. 

2.3 Matching in the moving- lattice representation 

The basic tool for a rigorous analysis of the form of the solutions of our matching 
constraints is the moving-lattice method of ref. [13] as reviewed in Appendix A 
below. Skipping the majority of details let us only note that for an analysis of this 
type, one of the recommended techniques seems to be the reduction of the problem 
to a — r plane. Preserving the definition of r = t(N, t) of Appendix A and replacing 
the definition of o = a(N, t) by another formula, 



<j = a(N,t,K,r) = n x 



[N + t] 2 + 



L — " 



(K + r) 



2£ 

we eliminate the coordinate R. A shortcoming of this approach is that our matching 
condition (12) transferred into the a — r plane has to be understood as the following 
quadratic equation for r, 

$i T 2 + 0J K ,r,t T + n K ,rA°) = (15) 
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where we abbreviated 



o" 2 cosh 2 <j H a sinh 2 cr 



This defines r = r K ^ t (N) on the lattice, the "motion" of which will be controlled 
not only by t and r but also, not so strongly, by K. Technically, the price to be 
paid is still reasonable - we get the closed form of the matching-compatible function 
r = r(er) as the two well known root formulae from eq. (15). Nevertheless, significant 
simplifications of the resulting picture may be mediated by the direct inspection of 
the equations in question. 



3 Solutions 

3.1 Matching equations in the a — t plane 

Building far-reaching analogies with the £ = special case would be misleading 
because the form of our matching constraint (12) is discontinuous in the limit £ — > 0. 
Thus, let us assume that £ ^ and study eq. (12) in its full-fledged form. Firstly, 
we abbreviate M. (a, r) = a sinh 2 a + r sin 2 r and Af(a, r) = o 2 cosh 2 a + r 2 cos 2 r 
and re-write our matching constraint (12) as the secular equation 

P <7, r, i?) = Q <t, r) + = 0, Q(a, r) = ) ! . (16) 

This enables us to formulate several obvious observations. 

[01] The shape of both the functions M.(o, r) and N(a, r) of two variables is easily 
deduced using their separability, X(a, r) = X(a, 0) + X(0, r), X — A4,Af. 

[02] The smoothness of the a— and r— dependence of the denominator J\f(a, r) 
facilitates also the determination of the shape of J 7 (a, r) = 1/Af(a, r). 
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[03] In a — r plane we may visualize the shape of the second fraction in (16) as a 
function which is constant along hyperbolas R(a, r) = \Jr 2 — a 2 = fixed. 

All these innocent-looking observations have several far-reaching though not always 
obvious consequences and form in fact a background for a rigorous analysis of the 
spectrum. 

3.2 A rigorous graphical interpretation of Q(<j,t) 

In more detail, observation [01] means that the surfaces defined by the two non- 
negative function(s) X(<j, 0) > have the form of the two only slightly different 
parabolic valleys with the same degenerate minimum (= zero) which coincides with 
the axis a = 0. The pertaining second components X(0,t) differ more from each 
other but both are adding a structurally similar perpendicular set of infinitely many 
parallel hills and valleys possessing a steadily increasing (though always finite) am- 
plitude. As an obvious result of the superposition, both the resulting surfaces X(a, t) 
cross the zero plane merely along certain ovals O* , and both of them only get neg- 
ative in their interior. 

The precise shape of these ovals (numbered by n — 0, 1, . . .) may fully rigorously 
be determined using the moving-lattice method (cf. Appendix A) but even without 
any use of the moving lattices the qualitative character of their shape is obvious 
and we may conclude that the zero lines of M. (a, r) and J\f(a, r) form the families 
of ovals O^ 4 and located within the stripes of r £ \(n + l/2)n, (n + l)ir] and 
t E [(n+ l/4)7r, (n + 3/4)7r], respectively. All of them are symmetric with respect to 
the reflection a — > — o and their size in the a direction increases with r. 

Examples of these structures may be found in both refs. [9] and [13] and another 
illustration appears in Figure 1 here. In fact, the Figure displays another surface 
Q(<j,t) = M.(o,t)/N{o,t) (within a narrow window of < Q < 0.05) but the 
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shape of the curve where M. vanishes (O^ 4 = V±) appears there clearly since the de- 
nominator J- {a, t) = 1/J\f(a, t) has its zeros, generically, elsewhere (cf. observation 
[02]). Besides the oval V 1 (and a part of Ofi 4 = V ) the picture displays another oval 
Of = D 1 of the zeros of the denominator M . Incidentally it lies within the chosen 
interval of r G (3, 7) and remains visible due to a numerical artifact of a spurious 
projection of an infinite discontinuity of the function T{a, r). 

Although the visibility of the discontinuities reflects just an imperfection of the 
graphical representation of the surface, in will prove useful in what follows. 

3.3 The role of the second component of V[a,r, R(a,r)} 

The presence of the subsurface generated by the second, R— dependent component 
T>(fy[R(a, r)] in eq. (16) does not violate the separation between a and r too much 
(cf. observation [03]). At the smallest absolute values of a we may safely return to 
the approximation of T>^ [i?(a, r)] by a function of a single variable, [tan 2 g(R)]/R fa 
[tan2£r/(L — £)}/t. This picture only becomes deformed, at the larger a, by being 
bent to the right, i.e., along hyperbolas R(<t,t) = constant. 

A clear understanding of the r— dependence of the whole surface V[a, r, R(<r, r)] 
will be obtained when we distinguish between the domain of the "small r" {where 
[tan2£r/(L — £)]/t ~ 2£/(L — £) is positive and virtually constant}, "medium r" 
{with the repeated quick growth of the curve [tan2£r/(L — £)]/r from minus infinity 
up to plus infinity within each interval of the constant length At = it (L — £)/2£} 
and "large r" {where the values of V^[R{a, r)] ~ 1/r become very small up to the 
very thin layers near the singularity hyperbolas H n }. Due to the local dominance 
of the latter singularities H n at any n — 0, 1, ... it is easy to imagine that the sign 
of the whole function V[a,r, R(a,r)] is positive and negative in their left and right 
vicinity, respectively. This "rule of thumb" enables us to deduce the sign of the whole 
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function T>[a, r, R(a, r)] in all our Figures. 

3.4 The left-moving hyperbolic discontinuities H n 

In the domain of the small shifts f < 1 the numerical values of the R— dependent 
component V^[R(a,r)] of eq. (16) remain almost constant and small. In this 
regime the above-mentioned "small-r" constraint r <C (L — £)/£ is not particularly 
restrictive so that the matching-compatible roots of equation V = remain very 
similar to their £ = predecessors in quite a large leftmost portion of the a — t 
plane. In our notation, the first few ovals Of L = V n of the zeros of the secular 
determinant stay only perturbatively shifted and deformed by an increase of I <C 1. 

With the growth of £ or A = £/(L — £) the leftmost discontinuity-hyperbola H of 
the surface V[a, r, R(a, r)] moves to the left and emerges in the right half of Figure 2 
where we choose the scale- independent parameter A = 11/40 which corresponds to 
£ = 11 L/51. This means that we are just leaving the domain of the small shifts £ <C 1 
so that the deformation of the nodal oval Of = V\ becomes perceivable, caused by 
the closeness of H to the £— independent discontinuity oval Di = inherited from 
the never-vanishing factor J 7 (a, r) = l/jV(cr, r). 

In a way which generalizes the illustrative Figure 2, each hyperbolic singularity 
H k (defined by the equation R(cr, t) — (L — £){k + 1/2)tt/£ with k — 0, 1, . . .) moves 
to the left with the growth of £ and A. Once it gets close to the N— th singularity oval 
D N _i, it touches it at a point with the coordinates cr^V^ = and = (N — 1/A)tt 

at the critical value A = 21/ (L - I) = (4k + 2) /{AN - 1) = \ { ^ ] of the shift. 

With the further growth of A the intersection of the hyperbola with the standing 
oval moves to the left and disappears, curiously enough, at a certain pair of points 
with the "last-contact" \a\ = > and r = < (N - 3/4)tt. The latter 

value lies slightly below the oval's end. Let us skip here the proof of this subtlety as 
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not too relevant. 

3.5 A completion of the list of the nodal lines 

We are now prepared to detect all the nodal curves of V[a, r, R(a, r)] and to de- 
termine their qualitative £— dependence in all the interval of £ G (0, L) and/or of 
A = X(£) G (0, oo). For the first inspiration we return to Figure 2 where the oval of 
zeros Of = V\ cannot be interpreted as a mere small perturbation of Of in spite 
of the fact that the singularity hyperbola Hq still did not touch the singularity oval 
Of = D 1 since A = 0.275 < = 2 /7 « 0.286. 

Still, the much more important observation made in Figure 2 concerns the emer- 
gence of the new curve W of the new zeros of the function V. At the chosen A 
this curve just entered Figure 2 at its right side. Our next Figure 3 confirms that 
the new nodal curve W moves to the left and gets deformed in a way reflecting the 
presence of a steep oval dip in Q(a, r) below r = 2n. We choose A = 0.355 which is 
still safely smaller than the lower estimate (4/c — 2)/ (4iV — 3) = 0.4 of the singularity 

(2. 1) 

hyperbola's "jumped-over" parameter X]2ut) ~ 0.403. 

The "next-step snapshot" of Figure 4 at A = 0.395 shows how the same dip 
deforms the shape of the oval Of = V 1 in the domain where the function of R is 
small. In the subsequent Figure 5 we finally see how the two curves of the zeros merge 

(2. 1) 

while a topologically new situation is created and sampled at A = 0.415 > \(2ut)- 

We may summarize that for the growing A the motion of the singular component 
tan2g(R) / R of our secular determinant T>(a, r) to the left gives a clear guide how 
to keep the £— dependence of its zero lines under full control. The emergence and 
the asymptotically hyperbolic shape of the new (and, in fact, not quite expected) 
non-oval curves W m of zeros follows immediately from the asymptotic smallness of 
the positive component Q(a,r) ~ 1/a 2 of V(a,r) at the larger \a\ 3> 1. 
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Due to the reasonably elementary character of the function T>(cr, r) we are able 
to understand that the pattern sampled by the Figures 2 - 5 is entirely universal. 
Always, step by step, the nodal ovals V n as well as their asymptotically hyperbolic 
nodal-line partners W m become deformed by the existence of the dip in the numerator 
function M(cr, r). 

Of course, after the hyperbola of singularities as well as its strongly deformed 
trailing nodal curve Wk "creep" over the fixed singularity oval Dj (as well as over 
its attached and strongly deformed zero curve Vj), the smoother shapes of both the 
nodal curves Wk and Vj are more or less recovered, and only their ordering remains 
permanently reversed. In spite of the apparent nonlinearity of the "creeping-over" 
effects, their details might again be analyzed algebraically, using an adapted version 
of the moving-lattice method of section 2.3. 

The most important reward compensating an increase in complexity of the latter 
recipe is that one becomes able to treat one of the two roots of eq. (15), say, as a 
"non-perturbative" solution at the small £. The most important example of its role 
are the hyperbolic nodal curves W m which move to the right in r with the decrease 
of £ and which disappear in infinity in the NSW limit of £ — > 0. 

4 Energies 

4.1 Graphical representation and classification 

On the background of the preceding material, what remains for us to do is a combi- 
nation of the above- described knowledge of the nodal lines of V[cr, r, R(a, r)] with the 
coupling-dependence constraint a x t = Z = (L — £) 2 g /2. A sample of the intersec- 
tions of this type (i.e., of a typical final solution) is offered in Figure 6 where A = 2.40 
is neither small nor large and where we choose Z = = 1.00 and Z = Z<® = 2.24 
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(= the critical "exceptional-point" value of ref. [9]) for illustration. The conclusions 
which are illustrated by this graph have a general validity: 



• We always have r > a > which means that all the real bound-state energies 
E n remain positive at Z > 0. 

• Some of the energies remain real at any value of Z > 0. They correspond to 
the intersections of the hyperbola a = Zjr with the hyperbolic nodal lines W m 
and may be called "stable", E = E$. 

• All the other energies E = E^f> correspond to the intersections of the hyper- 
bola a = Zjr with the nodal ovals At a sufficiently small Z the latter 
intersections remain real (see the line (a) with Z = 1 in Figure 6). 

• We may call the latter energies "unstable" as they merge in pairs and form 
complex conjugate doublets [4] beyond certain "exceptional-point" [24] values 
of £ and Z (illustration: the line (b) in Figure 6). 

The decomposition of the spectrum into its stable and unstable parts varies with £ or 
\ — (L — £) in an obvious manner. Hence, the stability pattern in the spectrum will 
be entirely different at the small and large A since in the former case the hyperbolic 
curves W n only generate the high-lying energies and vice versa. 

4.2 Numerical construction 

After all our previous detailed analysis of the qualitative features of the spectrum 
the numerical determination of energies becomes fully routine. Indeed, as long as we 
know t = Z /a, the rule r 2 — a 2 = R 2 leads immediately to the definition of 




(17) 
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In parallel to such an introduction of the closed function a = a(R) of R we may return 
once more to the recipe r = Z/a(R) and re-read it as another explicit definition of 
the second auxiliary function t(R) = Z/a(R) of R. 

In such a setting, the purely numerical determination of the bound-state energies 
is reduced to the search for the roots R n of eq. (12), i.e., of the zeros of the secular 
determinant 



V(R) = a 2 (R) cosh2a(R) + T 2 (R) cos 2 r(R) 



sin2 XR + 



+R [a(R) sinh 2 a(R) + t{R) sin 2 r{R)\ cos 2 A R (18) 

converted now in the function of the single variable R ~ \[E. An illustration of 
such a search is given in Figure 7 at a fixed choice of Z — 2. The quadruplet of 
the graphs of the secular determinant T>{R) = T>[a(R),r(R), R] is presented there at 
the four different values 1.25, 1.35, 1.45 and 1.55 of A (indicated along the vertical 
axis). In each of these graphs we magnified the vertical units near V(R) ~ and 
compressed them to a single point representing all the bigger values of \D{R) \ > s. In 
this way the picture samples the left-hand side of eq. (18) solely near its zeros. Our 
magnification of the vertical dimension marks these zeros by the virtually straight 
parts of the curve which are seen as practically perpendicular to the horizontal axis. 

The set of graphs in Figure 7 illustrates the A— dependence of the bound-state 
roots R n . We see that a pair of the unstable energies may merge and cease to be real 
after a fine-tuned growth of A. This illustrates the complexification of the unstable 
energies which is not caused by the growth of Z but rather by the growth of A. At 
the first sight this phenomenon looks like a paradox because we are now weakening 
the non-Hermiticity in fact. Fortunately, this paradox is still easily understood once 
we imagine (and check, say, in the spirit of Figures 2 or 3) that the growth of A 
"pushes" all the zeros (including of course also the nodal oval in question) to the 
left. Of course, this oval cannot get prolonged in the a direction because the function 
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M. (er, r) itself grows too quickly with a. This implies that the two real intersections 
of the oval with the hyperbola a = Z/r disappear because the latter curve grows to 
the left. 

In the light of an additional scaling in eq. (14) one may only admire the subtlety of 
the phenomenon, the verification of which very much profits from the exact solvability 
of the model. An independent confirmation of the absence of any contradictions may 
be also offered via a further simplification of mathematics. This inspires us to pay 
particular attention to the "most counterintuitive" limiting case where L — > oo. Such 
an analysis may be of an independent interest as it simulates, very roughly, the shape 
of the most popular antisymmetric and purely imaginary potential V(x) ~ ix 3 with 
real spectrum [22]. As long as this discussion already lies somewhat beyond the scope 
the present text, it is moved to the Appendix B. 

5 Conclusions 

After more than ten years of an intensive research many people now seem to be- 
lieve that we now better understand the key problems related to the so called 
VT— symmetric as well as to many other similar non-Hermitian models or, in the 
more rigorous terminology, to all the models where the metric remains nontrivial, 
rj 7^ I [25]. By the way, not all the related results are new. For example, Scholz 
et al [26] (inspired, presumably, by a few earlier mathematical as well as physical 
publications) studied the similar rj ^ I models more than ten years ago (!) and 
coined the name "quasi-Hermitian" for them. 

Still, one cannot deny that during the last cca seven years, a new and intensive 
excitement has been caused by the discoveries of the reality of the spectra in many 
VT— symmetric models. The emphasis of the research has been shifted, typically, 
to the explicit constructions of the charge C [20] or to the more detailed analysis of 
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what happens at the "exceptional" points where the reality of the spectrum is being 
lost [24, 27]. A few unusual features exhibited by our present model seem to offer 
another welcome and clear intuitive guidance in this area. 

We found our results interesting since the merger and subsequent spontaneous 
complexification of some "twin" pairs E^ ±twm ^ of the energies cannot be easily de- 
scribed within the usual textbook models where the metric is "trivial" , r)( trivia i) = I. 
It is also in this context where considerations based on our present model could 
lead to a deeper insight in the underlying mechanisms and mathematics, not only 
because our model is solvable but also because it proves able to provide different 
"twin-merging" patterns in the spectrum. Indeed, by the choice of the shape param- 
eter £ we may, up to a large extent, prescribe which particular excitations (say, in 
the low-lying spectrum) should remain robustly stable and which ones should form 
the unstable, fragile "twins" merging at some sufficiently large couplings gonial)- 

In the similar constructions and studies, one might feel hesitant whether his/her 
models should be simpler or more realistic. We believe that one should transfer the 
insight gained in the solvable models (like in the present one) to all the more realistic 
applications where just some approximate methods can be used. In this sense we 
already mentioned a parallelism between the role of the shift £ in our solvable model 
and of the exponent \i in the power-law potentials with VT — symmetry. 

It is encouraging to see that a certain nontrivial enrichment of the merging pat- 
tern has been detected, more or less in parallel, within the class of the power-law 
forces [27]. In this comparison, our present model's merit lies in its exact solvability. 
Definitely, it proves able to offer a comparably rich pattern of the mergers of the 
levels. 

This being said, the key phenomenological and "model-building" specific merit 
of our present new version of the VT — symmetric square-well model is still to be 
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seen in the "global" structure of its spectrum. There, one observes that the "fragile" 
and the "robust" levels seem to form the two sets which may be moved with respect 
to each other as a whole. Thus, the whole spectrum becomes "almost completely 
robust" in one extreme (which is "almost Hermitian") and "almost all fragile" in 
another extreme which is, near £ m 0, "maximally non-Hermitian" . 
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Figure captions 

Figure 1. A thin slice through the surface Q(<j,t) = M./M. 

Figure 2. A thin slice through the surface of the secular 
determinant V{a, r) at A = lj{L - £) = 0.275. 

Figure 3. Same as Figure 2, A = 0.355. 

Figure 4. Same as Figure 2, A = 0.395. 

Figure 5. Same as Figure 2, A = 0.415. 

Figure 6. Solutions at = 1.00 and = 2.24, intersections 
marked by circles, A = 2.4. 

Figure 7. Four re-scaled graphs of the function V(R). 

Figure 8. Graphical solution of eq. (28) (y = ujy/2, T = 1) 
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Appendix A: The method of moving lattice 

Secular eq. (12) and its descendants contain quickly oscillating trigonometric func- 
tions of arguments 2r and 2g. In the spirit of ref. [13] it makes sense to re-parametrize 
both these variables according to the rules 

T = r{N,t) = nN + nt, TV = 0,1,..., t e (0,1), 

g = ixK + nr, K = 0,1, . . . , r G (0, 1) 

which separate their "large" change (by an integer multiple of the period 2tt so that 
the trigonometric function itself remains unchanged) from a "small" change [within 
one period (0,2n)]. Thus, once we define 

\I/ = sin 2 r = sin 2 ir t, $ = cos 2 r = cos 2 n t , 

E = — tan 2 g — — tan n r , 

all our trigonometric functions in question become independent of both the integer 
variables. Thus, once we decide to work, say, in the a — R plane, we simply introduce 
a lattice £ t>r of points with coordinates 

a = <r(N,t) 



r(o-) vriV + Trt 
and 

where t G (0, 1) and r G (0, 1) are fixed while N — 0,1, . . . and K — 0,1, . . . remain 
variable. Our secular equation (12) then becomes more easily analyzed at the fixed 
t G (0, 1) and r G (0, 1) when it may be re-read as a simplified mapping a — > R with 

r, „ / x - ^ cosh 2 a + Z 2 $ t . . 

R = R t , r {<j) = ~ r X - ; — — — — i , 19 

a 6 sinn 2 a + a Z W t 
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i.e., Rt >r ~ S r \a\ at | a\ >> 1 while 



cr * 



at | er| <C 1 etc. In the subsequent step, remembering that the latter formulae hold 
on the lattice £(t, r) only, we must let this lattice move with the variation of t and/or 
r. Within each box numbered by the pair (N, K) of non-negative integers we would 
be able to re-derive all the qualitative geometric considerations of section 3 in an 
alternative, quantitative manner. 
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ip(x) 



Appendix B: Shallow well 

In the infinite-size limit L — > oo our model degenerates to a purely imaginary square 
well with asymptotic boundary conditions 

ip(±oo) = (20) 

and with the VT symmetric matching conditions in the origin, 

^(0) = 1, d x ^(0)=iG. (21) 

This means that we have the general solution 

cos k x + B sin k x, iG(Oi), k 2 = E, 

(22) 

(L + iN) exp(-ax), x E (£, oo), a 2 = iT 2 -k 2 , 
with T = y/g and with the the purely imaginary constant B = % G/k. 

B.l. Matching conditions at x = £ 

Let us split a = p + i q in its real and imaginary part with p, q > 0. This gives the 
rules p 2 + k 2 = q 2 and 2pq = T 2 , easily re-parameterized in terms of a single variable 

a, 

T 

p = qcosa, k = qsina, q = — ^=, a€.(0,£/2). (23) 

V 2 cos a 

The standard matching at the point of discontinuity is immediate, 

coski + B sinH = (L + iN) exp(-crf), 

-sink£ + B coskl = ~(L + iN) exp(-ai). 

k 

After we abbreviate a/k — —taxiQi, we get an elementary complex condition of the 
matching of logarithmic derivatives at x — £, 

G = -i k tan(fc + Q)l (24) 
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Its real part defines our first unknown parameter, G = G(a). Due to our normaliza- 
tion conventions, the imaginary part of the right-hand-side expression must vanish, 
Re[tan(fc + fi)£] = 0. An elementary re-arrangement of such an equation acquires the 
form of an elementary quadratic algebraic equation for X = tan k£. Its two explicit 
solutions read 



p + q 

IF' 



Xo 



p 



k 



or, after all the insertions, 



tan 



tan 



\Jl cosa( + ). 
£T sin 
72 



cos a 



(-) 



tan 



tan 



a 



(+) 



a 



(-) 



(25) 



(26) 



(27) 



These equations specify, in implicit manner, the two respective infinite series of the 
appropriately bounded real roots a = G (0, £/2). 



B.2. Energies 

For a G (0,£/2) the left-hand-side arguments in eqs. (26) and (27) run from zero to 
infinity and the functions oscillate infinitely many times from minus infinity to plus 
infinity. In contrast, the limited variation of the argument a makes both the right- 
hand side functions monotonic, very smooth and bounded, tan[(£ — a^)/2] G (1, oo) 
and tan[a( - )/2] £ (0, 1)- This indicates that our roots k = fe(a^) will all lie within 
well determined intervals, 

k ( n +) e (n+^,n+^j, n = 0,l,..., 

k^ G + -, m + 1^ m — 0,1, 

An additional merit of parametrization (23) lies in an unambiguous removal of the 
tangens operators from both eqs. (26) and (27). This gives 

4 +) = » + 5 - ^ *t> = m + 1 - ^, u « = f e(o,i). 
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After a change of notation with = u 2n and ) = u; 2n+ i, we may finally combine 
the latter two rules in the single secular equation 

2N + 2 - u N 



sin -uj n 



2 cos \-0J N ) N = 0,1,..., (28) 



In a graphical interpretation this equation represents an intersection of a tangens-like 
curve with the infinite family of parallel lines. This is illustrated in Figure 8. The 
equation generates, therefore, an infinite number of real roots ujn £ (0, 1) at all the 
non-negative integers N = 0, 1, . . .. The discrete spectrum is unbounded from above 
and remains constrained by the inequalities 

(i±m. < En < cAf+i)! (29) 

4 ~ ~ 4 

independently of the coupling T. 



B.3. Wave functions 

Equation (24) in combination with eqs. (26) and (27) determines the real parameter 

G = G (±) = — (30) 

q±p K J 

responsible for the behaviour of the wave functions near the origin [remember that 
B = iG/k in eq. (22)]. For its deeper analysis let us first introduce an auxiliary 
linear function of uj and N, 

2N + 2-u N (N+l/2 N+l\ 



^/R{u N ,N) = 



4T \ 2T ' 2T 



and re-interpret our secular eq. (28) as an algebraic quadratic equation with the 
unique positive solution, 

cos(^-uj n )= ] = . (31) 

V 2 / R(co N , N) + ^R\u N , N) + 1 
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This is an amended implicit definition of the sequence un. As long as the right hand 
side expression is very smooth and never exceeds one, the latter formula re-verifies 
that the root oj^ is always real and bounded as required. 

In the weak coupling regime (i.e., in the domain of the large and almost constant 
R 3> 1 with the small square- well height T or at the higher excitations), our new 
secular equation (31) gives a better picture of our bound-state parameters ojn = 
1 — 7] jv which all lie very close to one. The estimate 

£ 115 

-t]n — arcsm . — — — - — — + . . . 

2 IN R+ ^WTl 2R 48 R 3 

represents also a quickly convergent iterative algorithm for the efficient numerical 
evaluation of the roots u>n- One can conclude that in a way compatible with our a 
priori expectations, the value of p = p^ = Rea ~ q/2R is very close to zero and, as a 
consequence, the asymptotic decrease of our wave functions remains slow. We have 
q — qN — I mo " ~ k so that, asymptotically, our wave functions very much resemble 
free waves exp(— ikx). In the light of eq. (30) we have also ip(x) ~ exp(— ikx) near 
the origin. 

In the strong coupling regime (i.e., for very small R representing, say, the low- 
lying excitations in a deep well with T ^> 1) we get an alternative estimate 



£ . l 

- a; jv = arcsm \ / - 
4 V 2 



R - (VTTRZ - l)] w \ R - -R 2 + . . . < -. 



2 4 4 

In the extreme of R — > the present spectrum of energies moves towards (and 
precisely coincides with) the well known levels of the infinitely deep Hermitian square 
well of the same width / = (—£,£). In this sense, the "complex-rotation" transition 
from the Hermitian well to its present non-Hermitian VT symmetric alternative 
proves amazingly smooth. 

The wave functions exhibit the similar tendency. In the outer region, they are 
proportional to exp(— px) and decay very quickly since p = 0(R~ 1 ^ 2 ). The parameter 
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GrW becomes strongly superscript-dependent, 

G (+) = _Jf_ = (R 3 / 2 ), G<"> = -(q + p) = 0{R- 1 ' 2 ). 
q + p 

This means that in the interior domain of x G (—£,£), the wave functions with the 
superscript ^ and ^ become dominated by their spatially even and odd components 
cos kx and sin kx, respectively. In this sense, the superscript mimics (or at least keeps 
the trace of) the quantum number of the slightly broken spatial parity V . 



32 



